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ABSTRACT 

This paper considers some simple surface brightness (SB) estimates for galaxies in the 
Automated Plate Measuring Machine (APM) catalogue in order to derive homoge- 
neous SB data for a very large sample of faint galaxies. The isophotal magnitude and 
area are used to estimate the central surface brightness and total magnitude based 
on the assumption of an exponential SB profile. The surface brightness measurements 
are corrected for field effects on each UK Schmidt plate and the zero-point of each 
plate is adjusted to give a uniform sample of SB and total magnitude estimates over 
the whole survey. Results are obtained for 2.4 million galaxies with blue photographic 
magnitudes brighter than bj = 20.5 covering 4300 deg 2 in the region of the south 
galactic cap. Almost all galaxies in our sample have central surface brightness in the 
range 20 to 24 bj mag arcsec -2 . The SB measurements we obtain are compared to 
previous SB measurements and we find an acceptable level of error of ±0.2 bj mag 
arcsec -2 . The distribution of SB profiles is considered for different galaxy morpholo- 
gies for the bright APM galaxies. We find that early- type galaxies have more centrally 
concentrated profiles. 

Key words: surveys - galaxies: photometry - surface brightness. 



1 INTRODUCTION 

Surface brightness (SB) is one of the fundamental parame- 
ters describing a galaxy. It plays an important role in many 
diverse aspects of extragalactic astronomy and cosmology, 
from identifying the whole family of galaxy populations to 
modelling their different distributions and motions in the 
Universe. These aspects are generally related to the forma- 
tion and evolution of galaxies, as well as to the nature of 
large scale structures in the Universe. 

Unfortunately, accurate measurements of the SB of 
galaxies rely critically on high quality observations with per- 
fect sky conditions and adequate telescope time, so it is not 
easy to get high quality SB data for galaxies, especially 
for low surface brightness galaxies (LSBG). Large, homo- 
geneous data samples are even more difficult to construct, 
and so uniform SB data for galaxies is very scarce com- 
pared to other parameters, such as position, colour, redshift 
and morphology. Hitherto relatively few sets of SB measure- 
ments have been available; the total number of such galaxies 
is only a few thousand, and the samples have either concen- 
trated on LSBGs or been limited only to bright galaxies (e.g. 
de Jong & van der Kruit 1994; de Jong 1996; Jansen et al. 
2000; Heraudeau & Simien 1996; Lauberts & Valentijn 1989; 
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Impey et al. 1996; Morshidi-Esslinger et al. 1999, hereafter 
MDS). 

On the other hand it is not difHcuft to make a rough 
estimate of a gafaxy's SB using parameters refated to both 
the luminosity and size of a galaxy. These two kinds of basic 
parameters, such as magnitude and effective radius of galax- 
ies, are common in many catalogues. One can consequently 
obtain a larger sample of SB data of galaxies at the expense 
of some accuracy for individual objects. The statistical anal- 
ysis of a large number of galaxies will not sensitively depend 
on the accuracy of the SB measurement for each individual 
galaxy because reasonable errors in SB can be allowed for 
in the analysis. 

The APM catalogue includes both the magnitude and 
isophotal size for each galaxy so these can be used to esti- 
mate the SB of galaxies. Based on this strategy, some work 
has already used SB estimates from the APM survey. Cote 
et al. (1999) used the mean isophotal SB of APM galaxies 
to select the candidates for an HI survey of LSBGs. More 
recently Cross et al. (2001) calculated the effective SB of 
galaxies in the 2dF Galaxy Redshift Survey and estimated 
their bivariate brightness distribution, as well as the number 
and luminosity density of galaxies. 

In this paper, we consider a variety of SB estimates, and 
present a method for determining the central surface bright- 
ness, no, of galaxies in the APM Galaxy Survey. We apply a 
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field correction across each survey plate, and use a matching 
procedure over the whole survey area to yield a homogeneous 
SB sample of over 2 million galaxies with a wide magnitude 
range over a large volume of space. The results will be used 
in future papers to study the distribution of galaxies as a 
function of surface brightness, in particular their clustering 
properties. We describe the basic data in Section [| The 
principles behind the method are described in Section ^| 
together with an introduction of the raw APM parameters 
that relate to the SB of galaxies. The field correction and 
matching process are briefly outlined in Section [|, where we 
also discuss the uncertainty of our SB estimates and com- 
pare to other SB measurements. In Section |E| we consider 
the distributions of (j,o and the SB profiles of bright APM 
galaxies. Finally we summarize our results in Section pi 



2 THE APM CATALOGUE 

The APM Galaxy Survey is hitherto the largest galaxy sur- 
vey used for the investigation of the large-scale distribu- 
tion of galaxies. It contains over 3 million galaxy positions 
and apparent magnitudes, compiled by scanning 269 UK 
Schmidt Telescope (UKST) survey plates over the southern 
sky with the Automatic Plate Measuring (APM) machine. 
A subset of 2.4 million galaxies in 185 UKST fields south 
of declination —20°, compiled from photographic data of a 
consistently high quality and representing the first stage of 
the APM survey, was described by Maddox et al. (1990a,b). 
This subset covered the south galactic cap and is defined by 
the boundaries of UKST plates whose centres have galac- 
tic latitude between b = —30° and —40° (the exact limit 
depending on star number densities, and therefore on longi- 
tude), as well as declination. It is this subset of 185 plates, 
representing an area of 4300 square degrees, that we use in 
this paper. 

The APM identified images using a connected-pixel al- 
gorithm, selecting groups of 16 or more pixels having inten- 
sities above a fixed threshold optical density. With a mean 
pixel size of 0.51 arcsec (the plate scale varies slightly from 
plate to plate and across individual plates) , this corresponds 
to a minimum area of 4.3 arcsec 2 . The detection threshold 
t was set for each plate to be about twice the rms noise 
in the sky, leading to a mean threshold of 12 APM density 
units, corresponding to ~ 25 bj mag arcsec -2 above the sky 
background. The threshold varies within each plate due to 
vignetting and emulsion effects and also varies from plate to 
plate, because of the different photometric zero-points for 
each plate. The standard deviation in the threshold due to 
the variation in photometric zero points from plate to plate 
is 0.3 bj mag arcsec -2 . 

The images were classified as stars, galaxies, blended 
images and noise based on comparisons with the median 
stellar profile as a function of magnitude, as described in 
detail by Maddox et al. (1990a). In our present analysis we 
consider only images classified as galaxies. The total magni- 
tude was estimated from the isophotal magnitude and area 
of each detected image, assuming a Gaussian light profile. 
This process was optimized for faint images near the magni- 
tude limit where galaxy profiles are significantly affected by 
seeing. Different plates have different limiting magnitudes, 
and so to ensure a uniform limit, the final sample was fur- 



ther constrained by the magnitude limit of the poorest plate, 
20.5 bj. Rejecting images fainter than 20.5 bj in total mag- 
nitude produced the final APM Catalogue that we consider 
in this paper. 



3 MEASUREMENT OF SURFACE 
BRIGHTNESS 

3.1 APM data 

The APM Catalogue includes measurements of the isopho- 
tal flux and the image area within the detection isophote for 
each image. These parameters are used here to estimate the 
SB of each galaxy assuming an exponential light profile. The 
APM parameters are calculated using emulsion density D, 
rather than incident flux / measurements. Below the emul- 
sion saturation D m ax, it is fairly well approximated by a 
linear relationship, 

f af + fi for f<{D w -P)/a 

\ £> max for f>{D max -/3)/a { ' 

where a and P are constants for any particular part of the 
plate (Maddox et al. 1990a, b, Cawson et al. 1987), and the 
value of -D max is found to vary by 2% across the field (Mad- 
dox et al. 1990b). The incident flux /, is the sum of the sky 
flux / s ky, and the object flux / Q bj . So, the APM measures the 
sky-subtracted density -Dobj for particular pixels of objects, 
as 

D hi = D — -D s ky 

= af + 13 - (af sky + 0) (2) 
= a/obj 

where D s ky = a/sky + P is the local sky value which was 
obtained from a pre-scan and smoothing process (Maddox, 
1988). Obviously, the upper limit of D obi is, D sa t = O ma x — 
Dsky, for any particular position in a plate. The saturation 
density D sat corresponds typically to a magnitude surface 
brightness ~ 22.3 bj mag arcsec -2 above the sky, but this 
figure varies with the sky background level. Note that this 
saturation is quite severe: for a typical disk galaxy with a 
central SB of 21.4 bj mag arcsec -2 , the profile is saturated 
inside 0.8 disk scale- lengths. 

For each image, the APM measures the relative bright- 
ness by summing -Dobj for pixels above the threshold t ( ~ 25 
mag arcsec -2 ), 

/iso= D °*r ( 3 ) 

The threshold t is constant across the field of each plate, but 
varies from plate to plate according to the fluctuation of the 
sky density for individual plates (Maddox et al. 1990a,b). 
When working with APM data it is convenient to define an 
'isophotal APM magnitude', mi so , of an image as 

miso = + 2.5 log 10 I lso (4) 

This APM parameter is related to the actual isophotal mag- 
nitude, feiso, which is given by 

biso = Z - 2.51og 10 Jiso (5) 

where Z is the magnitude zero-point of the photographic 
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plate. Also, the isophotal size of the image is given by 

A iso = 1 ( 6 ) 

D obj >i 

/iso, t and Aiso are the three main APM parameters that we 
use to measure the SB profile of galaxies. In former work on 
APM galaxy survey, Maddox et al. (1990b) used these pa- 
rameters together with the assumption of Gaussian galaxy 
profiles to calculate a total APM magnitude m tot which 
was then calibrated to bj using CCD magnitudes. For faint 
galaxies near the magnitude limit of the APM survey this is 
a reasonable approximation, because the galaxy images are 
small and the profile shape is dominated by the seeing disc. 
But for brighter galaxies an exponential profile is a better 
approximation (see Section ^). In the following subsections, 
we mainly use these three fundamental APM parameters to 
investigate the SB of galaxies. 

Additionally, there are two other APM parameters that 
are related to the profile of galaxies. Both of them are density 
weighted parameters. One is a 2 , defined as 



T,D ohi >t( x2 + y 2 ) D °bj(x,y) 

-^iso 



(7) 



where, D OOJ (x, y) is the density at location x and y rela- 
tive to the image centroid. Theoretically, a 2 , together with 
Ai so can constrain the SB profile of the image, and judge 
whether it is best fit by a Gaussian, or exponential, or even 
an r 1 / 4 profile. Unfortunately, in practice this is not possi- 
ble, because a 2 changes by only about 10% between these 
three profiles, and this is the same order as the observa- 
tional error in a 2 (see Section ^). Additionally, a 2 can be 
seriously affected by other factors, such as the ellipticity 
and substructure of the image, and any other asymmetry of 
the profile which will tend to increase a 2 . This means that 
the measurement of a 2 is not as robust as Ai BO or Ji so , and 
so we do n ot u se a 2 in the measurement of SB in this paper. 

we consider the distribution of a 2 for bright 



3.2 



In Section 

galaxies which have a large image and hence a smaller er- 
ror in a 2 , and show that an exponential profile is actually a 
reasonable approximation to the SB profile of most galaxies. 

Another useful parameter is the ellipticity of the image, 
which is derived from density-weighted second moments. 
This is can be easily converted to the apparent axial ra- 
tio of the APM image (Maddox et al. 1990a; Lambas et al. 
1992), 

V 



' = I - ! • (8) 

apparent 

On the assumption that the areal profiles of an APM im- 
age have similar shapes, then e is independent of other pro- 
file parameters, and will not affect the measurement of SB. 
Howeve r, e will affect the estimation of a 2 as discussed in 
Section 5.2. 



3.2 Modelling the SB profile of galaxies 

If the image of a galaxy is represented by a circularly- 
symmetric pure exponential surface brightness profile, 



p(r) = po exp(— r/ro). 



(9) 



The central peak value po, and the characteristic size ro fully 
describe the image. In this case, ro is the exponential scale 
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Figure 1. Schematic representation of various quantities used 
in our estimation of /iq- The solid line represents the galaxy SB 
profile as a function of radius. D sa t is the maximum measured 
APM density, and r sa t is the radius where this density is reached. 
The shaded area above Z? sat represents the flux lost due to the 
saturation, /sat- The isophotal measurement threshold is shown 
by the horizontal dotted line, labeled t, and Ti so is the isophotal 
radius. The shaded area to the right of r- lso represents the flux lost 
due to the isophotal threshold, / ne id- The unshaded area below 
the solid line is the observed isophotal flux, 7i so . 



length. The total luminosity of the image, hot, is simply the 
integral of p(r). The integration can be divided into three 
parts which contribute to J t ot, 



•ftot — ^sat + -fiso + -ffi 



(10) 



where 7 sat , Ji so and 7fl e id represent the flux lost due to sat- 
uration, the observed flux and the unobserved flux beyond 
the isophotal size of the galaxy respectively. These can be 
written as: 



2tt 



p(r)rdr - 27rr sat D s , 



I iso = 2tt p(r)rdr + 27vr SSLt D B , 



and 



where 



f CO 

/ p(l 



/field = 27r / p(r)rdr, 



n so = V ^i so /7T = ro ln(poA) 



(11) 



(12) 



(13) 



(14) 



denotes the radius corresponding to the isophotal area for 
circular images and 



rsat 



ro ln(po/-D sa t) 



(15) 



denotes the radius corresponding to the saturated area for 
circular images. These quantities are illustrated schemati- 
cally in Figure Of these three parts, 7i so is the only directly 
observed quantity, and is given by J2d b >t ^* ob j ( ec l-|l)- The 
measured values of Iiso and Ai so can be used together with 
eq.(|l4|) and eq.([l5|) to find the parameters po and ro. The 
resulting equations are not soluble analytically, but are sim- 
ple to solve numerically. The actual procedure that we adopt 
is as follows: take an initial value of (pp, r' ); calculate the 
(^iso, IL>) from c q-©' eq. (|l^) and eq.|t|). Then compare 
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Figure 2. (a) Left: Simulated surface brightness estimates, fio (central SB), fi c (mean SB within r e ) and fi- lBO (mean SB within the 
detection isophote, eq. jlq)) for an isophote t = 25.4 mag arcsec -2 as a function of total magnitude, 6tot, for circularly-symmetric 
exponential profile galaxies. The value of ro is fixed in this example at 2.15 arcsec (a typical value for the scale length), and the threshold 
and saturation value are set as 12 and 200 APM density units per pixel above the sky level, (b) Right: The isophotal magnitude, 
T^iso, (dashed line) as a function of total magnitude, mtot, for circularly-symmetric exponential profile galaxies with ro = 2.15 arcsec, 
illustrating the effects of the isophotal threshold and saturation. The total magnitude relation is shown (as a solid line) for comparison. 
Results arc shown for both APM magnitude parameters (equations (W) and (19)) and for conventional magnitudes. 



with the observational APM parameters (Ai so , 7i S o), and cal- 
culate a feedback of (Ap' , Ar ). Several iterations will lead 
to the final stable result of (po, ro). The analysis can be 
adapted to non-circular images by assuming concentric el- 
liptical isophotes and using the APM ellipticity. 

The effects of seeing can also be included in this iterative 
process, but for most galaxies the correction is less than 0.1 
mag. Our estimates of po are based on 7i so and A[ so , and for 
images larger than about 100 pixels these measurements are 
not sensitive to the seeing. 



3.3 Surface Brightness 

The pair of parameters (po, ro) entirely describes the expo- 
nential profile of circularly-symmetric images, so they can 
be used to infer the value of any other characteristic SB pa- 
rameter for such a galaxy. In particular we can define the 
central SB in units of magnitudes per unit area, 

/io = Z-2.51og 10 (p ) (16) 



where Z is the zero-point for that particular UKST plate. 
We can also define an effective SB, 

He = fl + C (17) 

where the constant C depends on the exact definition of 
the effective SB. It can also be useful to define an effective 
radius, r c , which is the radius containing half the light of 
the galaxy. This is related to the scale size by r e = 1.678ro 
for a circularly-symmetric image with an exponential light 
profile. For the SB at r e , the constant is C = 1.822, and 
for the mean SB within r c , C = 1.124. Also the average SB 
within the isophotal area (for the APM scans of UKST blue 
plates this isophote is w 25 bj mag arcsec -2 ) can be written 
as 

^ iso = Z-2.51og 10 (J^j (18) 

where 7li so and 7i so are the corresponding isophotal size and 
brightness, derived directly from eq.(^|). Additionally, one 
can also calculate the 'total APM magnitude' from /tot, 

m to t = + 2.51og 10 (7 t ot) = + 2.51og 10 (27rporo), 

(19) 
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for circular images. This is related to the conventional total 
magnitude of the image, which is given by 



APM total magnitude parameter (m t ,) 
10 11 12 13 14 



Z-2.51og 10 (J t . 



Z - 2.51og 10 (27rp r ) 



(20) 



where Z is the magnitude zero-point. 

The relationship between these three SB measures can 
be seen in Figure |i[ In this Figure, we fix the value of ro to 
illustrate the dependencies between these SB parameters, 
giving ro a value of 2.15 arcsec, chosen to be fairly typical 
of the APM sample. Though the different SB measures are 
all derived from the same set of (po, ro) data, we see in 
Figure ^a that some galaxies can show much larger differ- 
ences between /j,i so and either fj,o or [i c than other galaxies. 
Figure ^b presents the isophotal magnitude, labelled mi ao , 
as a function of the total magnitude for the same value of 
ro- This is shown for both the standard magnitude and for 
the APM magnitude parameters of equations (^) and (|l9|). 
mi so can be significantly fainter than m to t, at both high 
and low SB. At high SB (bright total magnitude), the flux 
/sat is lost because of the emulsion saturation; at low SB 
(faint total magnitude) the flux 7fi c id is missed because of 
the isophotal limit. This also provides a rough estimate of 
how much no will be underestimated if these two parts of 
missing flux are not accounted for. The effect of saturation 
and isophotal threshold on the magnitud e is shown explicitly 
in Figure ^(a). As mentioned in Section |3.l| bright galaxies 
may be saturated over a large fraction of the image. This 
can be seen quantitatively in Figure ^|(b) , which shows the 
ratio of isophotal radius, ri so , and saturation radius, r sat , to 
the scale length, ro as a function of total magnitude. 

As described above, this approach only provides a rough 
measurement of the SB of galaxies. It does not use the real 
profile of different types of galaxies, which may be different 
from an exponential disk. It also neglects any internal struc- 
tures, such as arms, bars or a central bulge. Also, although 
we have rejected images that are most likely to be merged 
pairs, the automated image classification is not perfect, and 
there is a residual contamination of about 5-10% of merged 
images, whose profiles will not be well represented by our 
simple model. 

Although for any individual galaxy our SB measure- 
ment is unlikely to be very accurate, it is helpful to con- 
strain the shape of any given profile and it does allow us 
to make general comparisons of galaxy SB within the whole 
sample. Theoretically, this approach is suitable for any other 
profile which is specified by only two parameters, such as a 
Gaussian profile, 



p(r) = po exp (-(r/r ) 2 ) 



or an r 1//4 law profile, 



p(r) =p exp (-(r/r ) 



\l/4 



(21) 



(22) 



Without any other effective observational constraints on the 
profile shape an exponential profile is the most reasonable 
choice since it is a good representation for a majority of 
galaxies. As discussed in Section 5.2, we find that <r 2 as de- 
fined in eq.(^) is not very effective in distinguishing between 
different profiles. Additionally, the exponential profile has 
been a popular choice in earlier work, so it allows us to com- 
pare with other results. So, in this paper, we continue our 
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Figure 3. (a) Top: The change in total magnitude due to the 
isophotal threshold, Aiiimj, and saturation, Am S at, as a function 
of total magnitude, (b) Bottom: The ratio of isophotal radius, r; so , 
and saturation radius, r sa t, to the scale length, ro as a function 
of total magnitude. The value of ro is fixed at 2.15 arcsec, as for 
Figure 2. 



analysis and discussions under the assumption that galaxies 
have an exponential SB profile. 

Note that these estimates are in terms of the "APM" 
SB, which is based on APM density measurements. They 
are directly related to true SB only if there is a linear rela- 
tionship between the real flux and the APM density (eq.Q). 
Clearly the emulsion saturation violates this assumption, 
and CCD observations show that it is not exactly linear be- 
low saturation (Cawson et al. 1987). Our fitting technique 
approximates the saturation as a simple maximum APM 
density, and ignores any non-linearity below the saturation 
density. This non-linearity will not strongly influence the 
ordering of galaxies as a function of their SB, but more ac- 
curate measurements would have to allow for this factor. 
Furthermore, if the intrinsic SB is required for a particular 
analysis instead of the apparent SB, then the K-correction, 
cosmological dimming, and the correction of the apparent 
ellipticity of disk galaxies should also be taken into account. 



4 FIELD CORRECTION AND OVERALL 
PLATE MATCHING 

To provide a uniform sample of SB measurements over the 
whole APM galaxy survey area, we make two further sys- 
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Figure 5. Comparisons of the bj central surface brightne ss o f galaxies within the overlapping region of UKST fields F409 and F471 
(a) before any correction; (b) after field correction (Section 4.1); (c) afte r fur ther calibration of the magnitude zero-points by matching 
results for galaxies in overlapping regions of UKST plates (as in Section 4.2). The rms scatter in the differences are in turn 0.27, 0.25, 
0.25 mag arcsec -2 . The dashed line in each panel shows the line /xq(409) = /^o(471). 
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Figure 4. A contour plot of the field correction values, aver- 
aged over 185 plates. The dotted line shows A/xi = —0.015 mag 
arcsec -2 , the dash-dot line 0.0 mag, the dashed line 0.1 mag, the 
solid line 0.2 mag and the thick solid line 0.3 mag 

tematic corrections. First, we correct for the field response 
function to make the different parts of each UKST plate con- 
sistent from the central region to the edge. Then, we make 
an overall photometric zero-point adjustment to each plate, 
to give a uniform sample over the whole survey area. 

4.1 Field response corrections 

The variation in the response function across the field of each 
plate is caused by two main factors. One is the geometrical 
vignetting of the optical system of the telescope (UKST), 
that reduces the effective transmission off its axis. A theo- 
retical analysis (Dawe & Metcalfe 1982) shows the transmis- 



sion is very flat for images closer than 2.5° from the center 
of a plate, but quickly falls to ~ 80% of the central value at 
4° off axis. This reduction is same for each plate. Another 
loss in sensitivity is caused by the differences in the chemical 
desensitization at different parts of the UKST survey plates. 
For most of the survey plates, the original plate holder has 
an air gap between the plate and the filter, and the greater 
air circulation near the plate edges caused greater desen- 
sitization of the hypersensitized Ill-aJ emulsion. This effect 
causes a multiplicative reduction in the measured flux, which 
is quite similar to the geometrical vignetting, but the pattern 
actually varies for different plates, because the local weather 
conditions during individual observations were significantly 
different. 

Both of these effects can be considered as a change in 
the relation between APM densities and true flux. If two 
objects have same flux, /obj, but are located at different po- 
sitions on a plate, then the geometrical vignetting and the 
different desensitization lead to different observed APM den- 
sities -Dobj via eq.(|^). The constant j3 is essentially the fog 
level of the emulsion, which is constant, so the loss in sen- 
sitivity means that the slope a is not a constant across the 
plate. It is a function of position in the field, being reduced 
by the geometrical vignetting and emulsion desensitization 
from centre to the edge. Thus, both these effects can be 
considered as variations in a. 

Fortunately, a is only a function of position in the field, 
and does not depend on the specific properties of an image. 
Therefore, the correction process of Maddox et al. (1990b), 
which was devised to correct the total magnitude of images, 
can be simply applied here to correct the SB of galaxies. 
According to eq. (8) of Maddox et al. (1990b), the correction 
to the total APM magnitude is 

Ami = 2.51og 10 f^-J (23) 

where the subscript denotes the centre of a plate and the 
subscript 1 denotes the particular position of an object. Ac- 
cording to the relationship between the central SB and the 
total APM magnitude of eqs. (|l6|) and (H), we can see that 
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Ho will need the same correction as mtot: 



A/Lti = -2.51og 10 — 



«i 



(24) 



The value of A/ii is determined individually for each galaxy 
by estimating the APM sky density in the surrounding 10 arc 
minutes. The apparent variation in sky density is converted 
to an effective change in a using an empirical relation which 
was chosen to ensure that the number density of galaxies is 
uniform when averaged over many UKS plates. This correc- 
tion will be generally applicable for other SB measurements, 
except for the mean isophotal surface brightness, fiiso- Since 
the threshold t is fixed in terms of observed APM density, 
the actual threshold isophote varies over each field and a 
more complex correction would be required to ensure that 
/iiso is uniform. 

Most survey galaxies are located within the central 3° 
radius of a UKST plate, and so the correction is very small, 
typically < 0.1 magnitude. For galaxies outside this region 
the correction can increases up to 0.4 mag. See figure W for 
the pattern of averaged A/ii of field correction. See the fig- 
ures of Maddox et. al (1990b) for further details. 

Figures ^(a) to (b) show that the field correction slightly 
reduces the rms difference between the SB measurements for 
those galaxies imaged on two independent plates. F409 and 
F471 are a pair of contiguous plates with quite a large over- 
lapping area, and there is a total of 7960 identified galaxies 
with measurements from both the plates. For quite a few 
of these galaxies, the distance from the centre of one of the 
plates can be large, which means that uncorrected values 
of A/io can be quite large. Figure ^ (a) plots the before 
field correction and (b) after it. The rms scatter in these two 
plots is reduced from 0.27 to 0.25 mag. This improvement 
suggests that the rms A/j,q from the field correction is ~ 0.1. 



4.2 Matching calibrations in plate overlaps over 
the whole survey 

Variations in the emulsion sensitivity and observing condi- 
tions mean that each UKST survey plate has a different 
photometric zero-point. To produce a uniform set of galaxy 
SB estimates, we need to adjust the zero-point of each plate 
to force a match between galaxies in overlapping regions. In 
the APM survey paper II, Maddox et al. (1990b) described 
a scheme to match the APM total magnitudes over the sur- 
vey. In this paper we apply a similar procedure to match our 
SB estimates. 

The matching procedure uses the fact that each plate 
covers an area of sky that overlaps with its neighbours, so 
there are many galaxies that have independent measure- 
ments from two neighbouring plates. For a particular overlap 
the mean difference between the two measurements gives an 
estimate of the difference between the zero-points of the two 
plates. Taking all overlaps together we can find the set of 
plate zero-points which is most consistent with these plate- 
to-plate offsets. 

The main difference for this paper is that a galaxy SB 
profile is determined by two parameters (po, ro) as described 
in Section ^, and so both parameters must be matched. In 
principle, we could choose to match any two independent 
parameters derived from these two, but for simplicity we 
choose jj,o and & t ot- They are related to po and ro through 
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Figure 6. The frequency distribution of offsets in no for galax- 
ies in 809 overlapping fields. The dotted histogram shows the 
results before matching, with a = 0.43 mag arcsec -2 ; The solid 
line shows the distribution after matching, with a = 0.08 mag 
arcsec -2 . 
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Figure 7. The frequency distribution of offsets in mtot for galax- 
ies in 809 overlapping fields. The dotted histogram shows the re- 
sults before matching, with a = 0.31 mag; The solid line shows 
the distribution after matching, with a = 0.05 mag. 



eqs. (P) and @, while the zero-point in eq. ( p~6| ) represents 
the overall corrected uniform zero-point. Since we have as- 
sumed exponential profiles in this analysis, the magnitude 
zero-points are not the same as the Gaussian total magni- 
tudes that were applicable to galaxies near the plate limit, 
as used in previous work on the APM galaxy survey. 

Images were matched in the overlaps of every pair of 
plates with centres closer than 7°. For the 269 plates in the 
south galactic cap, this gives 809 distinct overlaps with an 
average of 1405 galaxies brighter than bj — 20.5 in each 
overlap. To allow for non-linearities which introduce a non- 
unit slope in the APM magnitude and SB measurements we 
used a linear fit to (fii — fij) vs (fit + fij)/^- Before matching, 
the variance in zero-point for the different UKST survey 
plates is very large, as can be seen in the offset of no for 
galaxies in the overlapping plate pairs. The distribution of 
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offsets jli — p,2 for all 809 overlaps, is shown in Figure y, and 
the rms scatter is 0.43 mag. After matching, this value is 
reduced to 0.08 mag. As shown in Figure]^, similar reduction 
in scatter is seen for m to t , where the rms before matching is 
0.31 mag, and after matching is 0.05 mag. 

For different plates, the zero-point corrections for 
and mtot are well correlated, but not exactly the same. The 
scatter about the correlation is caused by the different seeing 
on different plates, which leads to a different apparent ro- 
The matching process has effectively given ro for galaxies on 
each individual plate a systematic correction to a common 
"zero-point" . 

4.3 Error estimation 

For an individual galaxy, the deviation of the actual SB 
profile from an exponential profile is likely to be the main 
source of error in the SB estimate, but this effect is beyond 
the scope of this paper. If we assume a given profile, then 
we can discuss the "intrinsic" error in the APM SB caused 
by the observational errors in (iiso, ^4iso). 

The errors in (Ii so , Ai so ) come from a complex combi- 
nation of sources which we consider in two main categories. 
First is the observing process, which includes a mixture of 
weather conditions and instrumental effects, including the 
residual calibration differences between the survey plates 
that are not corrected by the field correction function and 
overall plate matching. This category also includes photon 
noise. Second are the errors from the APM scan itself. These 
errors can be estimated directly from repeat scans of se- 
lected plates. For 7i ao we find this source contributes only 
about 0.04 mag uncertainty, and for Ai so we find about 5%. 
These errors could be analytically transformed to fj,o follow- 
ing a standard propagation of errors procedure, but even if 
the resulting uncertainty in fj,o is slightly larger, it will be 
negligible compared to the observational errors. 

The overlapping areas of neighbouring plates provide a 
direct opportunity to estimate the combined errors in our 
SB estimates. Figure |B|(c) shows the residual errors for the 
overlapping field of F409 and F471, after applying the field 
correction and the overall matching correction, and it can 
be seen that there is still an rms scatter of about 0.25 mag 
between the two independent measurements for each galaxy. 
This implies the error in each fio will be ~ 0.18 mag. We 
have calculated the rms scatter for each one of the 809 plate 
overlaps, independently, and show the frequency distribution 
of these rms values in Figure |j| For most of these overlaps, 
the rms of Apo is less than 0.3 mag, and the mode is at 0.21 
mag, corresponding to a typical uncertainty of 0.15 mag for 
an individual SB measurement. 

If we consider all of the plate overlaps as a whole, the 
rms A/io over all 1,282,600 galaxy pairs is 0.23 mag, so the 
overall uncertainty for an individual SB measurement has 
an rms of 0.16. This value includes both the random er- 
ror on each SB measurement, and the rms residual differ- 
ence between the zero-points for all plates after applying 
the matching corrections. The difference between this rms, 
and the rms within each overlap suggests that the residual 
error in the zero-point matching has an rms of 0.07 mag per 
plate, consistent with the overlap residuals. 

Additionally, we have considered the error as a func- 
tion of fio- Table 1 lists the rms difference between the /xq 



o 




0.2 0.3 0.4 0.5 



rms of A/j, q (mag arcsec 2 ) 

Figure 8. Frequency distribution of scatters (rms) of A/in. of 
galaxy pairs for 809 overlapping fields. 



Table 1. The rms uncertainty in fio for galaxy pairs with different 
central surface brightnesses 



fio (mag arcsec 2 ) no. of galaxy pairs rms of A/io/^/2 



24.5 r- 


- 23.5 


40832 


0.22 


24.0 r- 


- 23.0 


202094 


0.19 


23.5 - 


- 22.5 


480355 


0.17 


23.0 - 


- 22.0 


639704 


0.16 


22.5 - 


- 21.5 


564329 


0.16 


22.0 - 


- 21.0 


373165 


0.16 


21.5 - 


- 20.5 


181962 


0.18 


21.0 - 


- 20.0 


62253 


0.21 


20.5 - 


- 19.5 


14258 


0.23 



measurements in overlaps for a series of fj,o ranges. Galaxies 
with no about 21.0 ~ 22.0 mag have the smallest error at 
0.16 mag. At both low and high surface brightness the error 
increases up to ~ 0.23. This tendency could also be seen in 
Figure ^(c). The reason for the increased error is that Ho de- 
pends on the profile of a galaxy (see Section [a.2| ), and when 
we estimate the total intensity of the galaxy, I to t (eqs. ( |lo[ ) 
to (|f])), only iiso is directly constrained by the observa- 
tional data, while /flow and I sa t, the other two contributions 
to /tot, involve extrapolations that are more extreme at low 
or high SB. 

So in summary, the overall "intrinsic" uncertainty in /in 
including measurement errors, and zero-point errors is typ- 
ically 0.2 mag. Comparing galaxies measured on one plate, 
the relative uncertainty is about 0.16 (also see the next sub- 
section) . 

4.4 Galaxy surface density and mean \i 

After applying the corrections to the individual plate data, 
we should be able to define galaxy catalogues selected with 
uniform SB criteria. We can make a simple test of the unifor- 
mity of the sample by calculating the mean SB for galaxies 
on each plate. Cosmic structure will introduce some intrin- 
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Figure 9. The distribution of mean fio for all plates in the survey. 
Circles show fields where the declination of the centre is < —20°, 
and crosses show fields where the declination of the centre is > 
-20°. 



sic variation in the mean SB in different directions on the 
sky, but this is likely to be small; indeed the variation in 
galaxy numbers from plate to plate at b.j — 20 is only 12% 
(Maddox, Efstathiou & Sutherland 1996). Figure ^| shows 
a plot of the mean fio for each plate in the survey, before 
and after applying the matching corrections. Using all fields, 
the rms scatter before correction is 0.25 mag, and after cor- 
rection it reduces to 0.095 mag. However, it is clear that 
fields centred with declination > —20° have a larger scatter, 
even after correction. This is partly due to the fact that the 
scanned plates with declination > —20° are a rather hetero- 
geneous mix of original plates and non-survey plate copies, 
and partly due to the significantly large galactic extinction 
in this part of sky (see Schlegel, Finkbeiner & Davis 1998). 
In principle it should be possible to correct these problems, 
but for simplicity we have simply excluded the fields with 
declination > —20° in the subsequent analysis. 

A further simple check on the uniformity is to compare 
the average number of galaxies per plate for SB selected 
samples. Figure [lo| shows the distribution of galaxies on the 
sky for three subsamples selected on SB from high to low SB 
respectively. The ranges in no were chosen to give samples 
with roughly the same number of galaxies in each subsample. 
Note that the high and low SB samples are selected from the 
tails of the SB distribution (see Section pM| ) and so are most 
sensitive to errors in the plate matching. 

For the 185 plates there are no plates which show a 
clearly discrepant number of galaxies, suggesting that our 
matching procedure does indeed produce a uniform set of 
SB measurements. This improvement can also be seen in 
figure |ll|, which plots the numbers of galaxies in each plate 
for these three subsamples before and after matching cor- 
rection. After matching, the fractional rms scatters in the 
number of galaxies on each plate are 0.27, 0.13, 0.17 re- 
spectively for the three samples. These are close to what is 
expected for randomly selected 5° x 5° patches of sky. 

The higher rms for the high SB sample reflects the 
higher clustering amplitude of high SB galaxies, which is 
clearly apparent in Figure [l(]. We plan to quantify these 
differences in in clustering in a future paper. 



(a) 



(b) 



(c) 






Figure 10. The distribution of SB selected subsamples of galax- 
ies projected on the sky: (a) 20.5 < /Ltrj < 21.2 bj mag arcsec -2 ; 
(b) 22.0 < no < 22.2 bj mag arcsec" 2 ; (c) 23.0 < < 24.5 bj 
mag arcsec -2 . Each panel shows a greyscale map shaded accord- 
ing to the surface density of galaxies in 14' X 14' cells over the 
SGP. White represents galaxies per cell, and black represents 
4 times the mean surface density. The dashed lines show right 
ascension and declination. 



4.5 Checking the APM results and photometric 
calibration with deeper CCD data 

It is important to test the accuracy of the APM surface pho- 
tometry, particularly the photometric calibration, through 
a comparison with well-calibrated higher-resolution CCD 
data. To this end, new CCD surface photometry has been 
performed on faint (bj ~ 18 to 22 mag) galaxies in the re- 
gion around the Hubble Deep Field South (HDF-S), which 
lies in the area of the available APM data (at galactic coor- 
dinates (/,&) = (328°, -49°)). Much well-calibrated public 
imaging data are available around the HDF-S, and ground- 
based HDF-S CCD observations cover a sufficient area to 
provide the numbers of galaxies required for these tests. The 
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Figure 11. The number of SB selected subsamples of galaxies 
per UKST survey plate's central 5° X 5° region, before and after 
matching correction. Their rms scatters, normalized by their av- 
erage numbers of galaxies, are 0.40, 0.20, 0.35 separately before 
matching and are reduced to 0.29, 0.13, 0.19 after it. 



Goddard Space Flight Center HDF-S public data (Teplitz 
et al. 2001; Palunas et al. 2000) were selected for the task 
because they combine a large image area, standard broad- 
band photometric passbands and good photometric calibra- 
tions. The Goddard HDF-S data were recorded with the Big 
Throughput Camera on the 4-metre Blanco telescope at the 
Cerro-Tololo Inter- American Observatory (Wittman et al. 
1998) and cover a 47 arcmin by 46 arcmin area. An inde- 
pendent measurement found that the seeing of the blue data 
was 1.7 arcsec in the central 15 arcmin by 15 arcmin region; 
although not good, this is comparable to or better than the 
photographic APM data. Pixels were 0.43 arcsec in width. 

We measured the SB profiles of galaxies in this region 
through the following steps. Firstly, images were detected in 



the B-band CCD data using the SExtractor package (Bertin 
1998; Bertin & Arnouts 1996). This provided image centroid 
positions, ellipticities and orientations. For each galaxy im- 
age, mean surface brightnesses were computed in a series 
of concentric elliptical annuli defined by these SExtractor 
image parameters. Then the surface brightness was plotted 
against semi-major axis for each image and an exponen- 
tial light profile fitted, defining the range in semi-major axis 
for these fits by visual inspection. Care was taken to avoid 
fitting to the core of the image where seeing effects were 
significant. This provided a central surface brightness and 
exponential scale length for each selected galaxy. Finally, 
the surface brightness results were transformed from the in- 
strumental blue magnitude system to the photographic bj 
system using the Goddard HDF-S transformation and the 
standard B to bj equation of Blair & Gilmore (1982). This 
used colours measured within identical elliptical apertures in 
the B and V data. The overall uncertainty in the individual 
central surface brightness measurements is about 0.1 mag. 

An independent check of these results was made us- 
ing other public ground-based HDF-S data. SExtractor to- 
tal magnitudes were computed using the Anglo- Australian 
Observatory public CCD observations for the region imme- 
diately around the WFPC2 field. These were transformed 
to the bj band and compared with the total magnitudes 
under the fitted exponential profiles derived from the God- 
dard HDF-S data. Six galaxies with Goddard photometry 
were identified in the AAO field. The difference in the total 
magnitudes was +0.032 ±0.030 mag, confirming the validity 
of the results obtained from the Goddard public data. The 
surface photometry from the Goddard data can therefore be 
used to check the APM photographic results. 

The HDF-S CCD surface photometry gave central sur- 
face brightnesses [io and scale lengths ro for a large sample 
of galaxies with bj — 17.5 to 22 for which the exponen- 
tial profile was judged to be an adequate fit. Of this sam- 
ple, 24 objects appeared in the APM Galaxy Catalogue. Of 
these 24, two galaxies were rejected from further consider- 
ation because their angular sizes were so small (isophotal 
areas smaller than 100 pixels) that seeing was likely to af- 
fect appreciably the results from the APM data. This left 22 
galaxies in the region around the HDF-S that had surface 
photometry both from the APM data and from the CCD 
observations. These galaxies were used to test the reliability 
of the surface brightness results derived from the APM data. 

Figure ^| shows the relationship between the CCD no 
and our APM measurements. The rms scatter is 0.29 mag 
arcsec -2 , and the difference in their zero-points is 0.26 mag 
arcsec -2 . This dispersion is larger than expected from sim- 
ply adding the uncertainty of 0.18 in the APM measure- 
ments in quadrature with 0.1 from the CCD measurements. 
The extra scatter is probably due to the fact that we have 
not corrected the non-linear correspondence between flux 
and the emulsion density of the survey plates. Note also that 
nearly all of the galaxies in this sample have fio < 22 mag 
arcsec -2 , which is brighter than the saturation level of the 
UKST emulsion and thus the APM measurements depend 
on extrapolated galaxy profiles. 

The scatter could be reduced by restricting the range 
in radius used to fit the CCD data so that it matches more 
closely the range available in the APM data. However, this 
would not be a true representation of the uncertainties in 
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Figure 12. Comparison of central surface brightness of HDF-S 
and this paper. The rms scatter is 0.29 mag arcsec -2 . 



the APM SB estimates. We want to compare the APM re- 
sults with results from CCD surface photometry that have 
been obtained using standard surface photometry methods. 
Note that there will also be systematic biases correlated with 
morphological type. Since we have assumed an exponential 
profile, a galaxy with an r 1 ^ 4 pro file will not be well repre- 
sented by our model (see Section 5.2). 

In summary, although the APM measurements have 
relatively large uncertainties, they are adequate to classify 
galaxies into broad surface brightness categories. 



4.6 Comparison with other SB measurements of 
UKST survey plates 

Morshidi-Esslinger et al. (1999, hereafter MDS) made a sur- 
vey of a substantial number of UKST survey plates in the 
south galactic cap in an attempt to identify LSBGs with 
large angular size. They used the APM to measure profiles of 
nearly 3000 LSBGs and estimated the central surface bright- 
ness by extrapolating the best fit exponential profile. Their 
SB estimates are based on the same original plates as our 
data, and they use the same scanning machine. They also 
make the same assumption of exponential galaxy profiles, 
so the MDS survey is not a completely independent data 
set. Nevertheless it provides a useful comparison for our es- 
timates. 

There are two main differences between our measure- 
ments and those of MDS. First, they fitted exponential pro- 
files to up to 8 isophotal areas Ao, Ai, Aj - see the details 
in Maddox et al. (1990a) and MDS (1999) - to solve for the 
(po, ro) parameters, whereas we used only Ao to calculate 
these values. The second difference is in the calibration pro- 
cedure. MDS's calibration of SB depends on zero-points for 
total magnitudes from other works for fields where this is 
available, and uses the magnitudes of stars in plate overlaps 
for fields with no external calibration. Our calibration de- 



pends mainly on matching fio directly for the galaxies plate 
overlaps, with an overall calibration from external photom- 
etry. 

Figure [l3| compares our estimates of no with 2185 MDS 
galaxies that can be identified in our sample. The unidenti- 
fied galaxies are mainly classified as 'merged' images in the 
APM galaxy survey, and thus not included in our catalogue. 
There is a very good correlation between the two estimates, 
but the scatter is rather large, with an rms of 0.35 mag. 
To investigate the cause of this dispersion, we have exam- 
ined data from three fields (F157, F358 and F482) having 
more than 65 MDS galaxies in each, and have plotted them 
with different symbols in Figure |l^. It is clear that the rms 
scatter for each plate is very small, only 0.16, 0.15 and 0.14 
respectively. These scatters are slightly smaller than the typ- 
ical rms uncertainties estimated from the APM plate over- 
laps. So, even if we assume that the MDS measurements are 
correct, and all of the scatter is due to the approximation 
involved in our simple method of estimating /xo, the errors 
introduced are smaller than the observational errors associ- 
ated with measuring different UKST plates. 

These three plates were deliberately chosen to lie at 
the extremes of the distribution plotted in Figure [H^, and 
it is clear that the zero-points of these three plates are very 
different. The zero-point for field F358, is almost 1 mag away 
from other two plates. If our zero-point was in error by 1 
mag, it would be clearly visible as a discrepant point in 
Figure []. Also the area covered by field F358 would show up 
as a large overdensity in Figure |lO|(a), and an underdensity 
in Figure [ji](c). None of these tests suggests there is any 
significant error in our zero-point for this plate, and so we 
believe the MDS zero-point is in error. Assuming that the 
rms err or o n our measurements is 0.16 mag as discussed in 
Section 4.3, this suggests that the MDS zero-points have an 



rms of 0.31 mag, and so F358 would be a 3 sigma deviation, 
which is not unexpected in a sample of around 100 plates. 

MDS estimate that the uncertainty in their calibration 
is about 0.1 mag for plates with external calibration, and 
0.2 for plates where they have used overlaps. There are two 
possible sources for the larger errors that we find. First, they 
do not apply any sort of vignetting or desensitization cor- 
rection. Since their zero-points rely on edge matching and 
small-area CCD calibration data, this will introduce an error 
in the plate zero-points with an rms of order 0.1-0.2 mag. 
Second, MDS use the magnitudes of faint stars in the plate 
overlaps to estimate the zero-points on plates with no ex- 
ternal calibration. Even faint stars tend to have a higher SB 
than galaxies, and this is likely to introduce a further bias. 
These two effects together might account for the observed 
0.31 mag rms difference. This comparison shows that cal- 
ibrating the zero-points of photographic data is a difficult 
problem that is still without a satisfactory solution. 

Impey et al. (1996, hereafter ISIB) scanned 24 UKST 
survey plates to search for LSBGs with fio > 22 mag 
arcsec -2 . Some of these plates are covered by the APM 
data used in this paper, and we find 202 galaxies in com- 
mon with their data. Plots comparing the our magnitudes 
and surface bri ghtn esses with the ISIB measurements are 
shown in Figure [15l There is a good correlation between the 
ISIB estimates and ours, but there are several differences in 
their analysis and definition of parameters which preclude a 
quantitative comparison with their data. 
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Figure 15. Comparison of total magnitude, central surface brightness and effective surface brightness of ISIB and this paper. Note that 
all data from ISIB are in the B band, while our APM results are in the bj photographic band. 
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Figure 13. A comparison of central surface brightnesses of MDS 
and this paper. The dispersion is cr = 0.35 mag arcsec -2 . To be 
consistent with MDS(1999), we have applied a correction to our 
fig to account for the ellipticities of galaxies. 
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Figure 14. The same as figure |13| but showing data from only 
three independent plates which have more than 65 MDS galaxies 
on each. The three fields are: F157, which has an rms scatter of 
0.16 mag arcsec -2 , is plotted with diagonal crosses; F358, which 
has an rms scatter of 0.15 mag arcsec -2 , is plotted with verti- 
cal/horizontal crosses; F482, which has an rms scatter of a = 0.14 
mag arcsec -2 , is plotted with open circles. 



5.1 The Distribution of Surface Brightnesses 

Figure [Hi] shows a histogram of the frequency distribution of 
fxo. The distribution is centred on = 22.2 mag arcsec -2 , 
and extends to fig ~ 20 mag arcsec -2 at the high SB end, 
and to fio — 24.5 mag arcsec -2 at the low SB end. The star- 
galaxy separation criteria (Maddox et al. 1990a) mean that 
some very high SB images will be rejected from the galaxy 
sample because their profile is so similar to stars. However, 
the rapid drop in the apparent number of high SB galaxies 
will be due to their small isophotal areas and a real drop 
in the space density of high SB galaxies (as shown by Cross 
et al. 2001). The lack of very low SB galaxies is less clearly 
a reflection of the true distribution of galaxy SBs, since a 
galaxy must have at least 4 square arcseconds brighter than 



25 mag arcsec -2 to be selected in the APM galaxy sam- 
ple. The distribution is further constrained by the 20.5 mag 
limit on the Gaussian profile magnitudes which will also bias 
against low SB galaxies. Cross et al. (2001) have used the 
2dF Galaxy Redshift Survey to estimate the distribution of 
intrinsic surface brightnesses, and conclude that, although 
the APM selection criteria do exclude low SB galaxies, the 
missing galaxies would not contribute significantly to the 
overall luminosity density. We will consider how our selec- 
tion criteria affect the observed distribution in a future pa- 
per. 

It is clear from Figure [iij] that the fio subsamples plot- 
ted in Figure |Io|(a) and [To](c) represent the galaxies at the 
extremes of SB in the APM sample. This means that the 
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Figure 16. The frequency distribution of the apparent SB, fio, 
for all galaxies in the survey. 

distribution of galaxies in the sub-samples is very sensitive 
to any residual errors in the uniformity of the /xo measure- 
ments, and so the apparent uniformity of the galaxies in 
t igurc 0(a) and 0(c) show that our corrections have given 
us a reliable and uniform set of measurements. 

5.2 The profiles of APM bright galaxies 

In this Section, we use galaxies from the APM Bright Galaxy 
Catalogue (Loveday 1996; hereafter APMBGC) to investi- 
gate how the profile of a galaxy depends on its morpholog- 
ical type, and test if the assumption of an exponential pro- 
file is valid. The galaxies in the APMBGC have bj < 16.44 
and have a relatively large angular size (riso) and large a 2 . 
Although <7 2 is not very sensitive to changes in the pro- 
file of galaxies, and the observational error is quite large, 
this parameter does contain useful information about the 
SB profile. Furthermore, the APMBGC galaxies have been 
morphologically classified by visual inspection of the UKS 
plates and we have analyzed the distribution of a 2 separately 
for both early and late type galaxies. 

We calculated the (po, ro) of each galaxy, assuming an 
exponential SB profile as in Section ^. Since we also know 
the apparent ellipticity e for each galaxy (from the density 
weighted second moments), and the values of the detection 
threshold and saturation of the plate, we can simulate the 
observed APM density for each pixel. Any other observa- 
tional parameters can then be estimated from this expo- 
nential model image. In particular, we have estimated the 
value of <t 2 , and compare this to the actual observed value. 
Obviously, if the original image is exactly an exponential 
profile, and there are no other observational errors, then 
the observed value, a Q h s will be exactly equal to that from 
the exponential model a cxp , so the ratio <r b s /o"oxp should 
be exactly unity. To test the sensitivity of a to variations 
in the SB profile, we created simulated images that follow 
the r 1//4 law profile (eq. |22| , and analysed them using the 
assumption of an exponential profile. In this case we found 
that o-obs/o" OX p — 0.9, with a small variation according to 
their real central surface brightness and apparent elliptic- 
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Figure 17. The frequency distribution of o" o bs/ 'exp values for 
APM bright galaxies. The dotted line shows E and SO galaxies; the 
solid line shows Spiral and Irr . An exponential profile would have 
CT obs/°"exp = 1-0, an r 1 / 4 law profile would have <T b s /cr eX p = 0.9, 
and a Gaussian profile would have cTobs/cXexp = 1.1. 

ity. The same analysis for Gaussian profiles gave a ratio of 

O" obs /Oexp — 1.1- 

In Figure [I?], we plot the frequency distribution of 
Oobs/oexp for early and late galaxies from the APMBGC. 
There is a large dispersion in <T b s /o" oxp , even for only late 
type galaxies. Nevertheless, the peak of the distribution is 
very close to 1.0, which implies that the exponential profile 
is a reasonable approximation for most galaxies (77% of the 
sample has 0.9 < o- bs/o" x P < 1.1). For early type galaxies 
(E and SO), the dispersion is larger than for late type galax- 
ies, and most of them have smaller (T bs/o"cxp values that are 
more consistent with an r 1//4 law profile. This shows that, al- 
though the APM measurements span a rather limited range 
of SB, the a parameter can distinguish between early and 
late-type galaxies. 

For fainter APM galaxies, the observational uncertainty 
on a 2 will be much larger, and the small differences between 
the r 1 / 4 , exponential and Gaussian assumptions mean that 
a 2 is not effective in constraining the SB profile. 



6 SUMMARY 

We have described a simple method to estimate the central 
surface brightness, fio, of galaxies, allowing for saturation 
and isophotal effects, and applied it to all the galaxies in the 
APM galaxy survey. We have corrected the /io estimates for 
the effects of vignetting and differential desensitization over 
the field of each UKS plate. We have also used the plate 
overlap areas to match the individual plate measurements, 
and so provide homogeneous fio estimates for all the APM 
galaxies. Internal and external checks of the uncertainties in 
/zo suggest that the uncertainty for an individual measure- 
ment is about 0.16 mag. 

For very bright APM galaxies, the second radial mo- 
ment of the SB profile, a 2 , shows that late type galaxies are 
consistent with an exponential profile, whereas early type 
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galaxies have a higher mean SB, consistent with an r 1//4 
law. 

This uniform set of SB estimates for the APM galaxies 
will be used in later papers to select samples of high and low 
SB galaxies and analyze differences between their distribu- 
tions. 
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